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We propose a quantum optomechanical implementation of the quantum mechanics free sub- 
system (QMFS) proposed by Tsang and Caves [Phys. Rev. Lett. 105, 123601 (2010); Phys. 
Rev. X 2, 031016 (2012)] based on the mechanical properties of a Bose- Einstein condensate trapped 
in an optical lattice. In that scheme the effective mass of the atomic condensate is modulated by an 
optical lattice and can become negative, resulting in a negative-frequency optomechanical oscillator, 
negative environment temperature, and optomechanical properties opposite to those of a positive 
mass system. We also discuss possible protocols to apply this scheme to the back-action evading 
measurement of feeble fields. 
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A fundamental issue faced by quantum metrology is 
the need to circumvent measurement back-action noise. 
Many of the early ideas were developed in the context 
of gravitational wave detection [l|-i4j , with groundbreak- 
ing advances such as back-action evading - or quantum 
non-demolition (QND) - measurements, as well as to the 
realization of the potential benefits of using non-classical 
fields, such as squeezed states, to locate the unavoidable 
quantum fluctuations where they would not significantly 
perturb the measurement. Recently these ideas have wit- 
nessed a significant new development by groups around 
C. M. Caves [5| and E. Polzik [^, who realized that it 
is sometimes possible to isolate quantum-mechanics free 
subsystems (QMFS) of a quantum system. All observ- 
ables in these subsystems are by construction QND ob- 
servables. As such they may find a number of applications 
in the detection of feeble forces and fields, including op- 
tomechanical sensing, magnetometry and gravitational 
wave detection. 

As discussed by Tsang and Caves [5| , a simple setup to 
implement a QMFS comprises two harmonic oscillators 
of identical frequencies and opposite masses, 
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Since [Q,n] = this means that the dynamical pair of 
observables {Q,!!} formed by the collective position Q 
and relative momentum 11 form a QMFS - and likewise 
for the pair {$,7^}. An implementation of this idea in- 
volving two spin ensembles oppositely polarized along an 



external magnetic field was recently demonstrated [y], 
and a system involving a cavity mode driven by a blue- 
detuned pump laser has also been advanced l7|. Here we 
propose a quantum optomechanical setup [8| where the 
effective negative mass system is produced by a Bose- 
Einstein condensate trapped in an optical lattice. 

Cavity optomechanical systems based on the collective 
motion of Bose-Einstein condensates (BEC) [9( or non- 
degenerate ultracold atomic gases [lO| have proven to be 
particularly well suited to demonstrate a number of quan- 
tum effects, including the observation of the quantum 
back-action of position measurements |10| , the asymme- 
try in the power spectrum of displacement noise due to 
the non-cummutating nature of boson creation and an- 



nihilation [11[, as well as the optomechanical cooling of 
a collective motional mode of an atomic ensemble down 
to the quantum regime [l2| . These experiments pave the 
way to promising ultracold atoms-based quantum metrol- 
ogy schemes. 

The concept of effective masses - which can in princi- 
ple be negative as well as positive - is familiar from solid- 
state physics where it has proven useful in describing the 
motion of electrons in non- ideal lattice potentials [13| . 
Not surprisingly, that idea has recently been expanded 
to describe aspects of the dynamics of ultracold atoms 
in optical lattices [IJ, |l5|. We show in this paper that 



when combined with BEC-based quantum optomechan- 
ics, it leads to the realization of QMFS that may find 
applications in the detection of feeble fields and forces. 

Consider for concreteness a scalar atomic BEC con- 
fined by both an optical lattice potential Vb cos^(fci,a;) 
of periodicity 27: /k^ and an external trapping potential 
U{x) that is taken to be slowly varying over the lattice 
period. Restricting the description to one dimension for 
simplicity, the condensate is described by the Hamilto- 
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where '^(x) is the bosonic field operator of the atomic sys- 



tern and we have neglected inter-atomic collisions. As- 
suming that the lattice potential is sufficiently shallow 
that we are far from the Mott insulator transition [Ig 
we proceed by expanding '^(x) in terms of a complete 
set of basis Bloch functions as 
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where n and q label the band index and the quasi- 
momentum, respectively, and d„^g arc the associated bo- 
son annihilation operators. In the following wc assume 
that the condensate is properly described by the product 
of an envelope that varies slowly over the period of the 
lattice, and is characterized by a central wave vector qo, 
and Bloch functions that capture the rapid oscillations of 
the condensate caused by optical lattice. We then have 
approximately [17[ 
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where the mode functions (f>n,qo i^) capture the density 
oscillations and 
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Here we have introduced the slowly varying bosonic op- 
erators (on the scale of the lattice period 2-K/k]^) 
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that describe the dynamics of the condensate envelope in 
the trapping potential U{x), with 
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Applying the effective mass method [l3| yields then 
the effective Hamiltonian describing the condensate in 
the envelope representation as 
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where the effect of the lattice potential is now included 
in the single-particle effective mass 



m* = hyCiqo). 
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Here ^^(go) and £"(.10) are the first and second-order 
derivatives of the nth-band Bloch energy with respect 
to the quasimomentum, evaluated at go- Figure [T] shows 
that the gradient term £,^ (go) vanishes at the center (go = 
0) an edges (go = ±Al) of the first Brillouin zone. 

The effective mass is negative at the zone edges for 
odd-ri bands, and at the zone center for even-n bands. 
Note also that the potential U{x) is shifted by the Bloch 
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FIG. 1. Top: Bloch energy £{qo) (a), its derivative (b), and 
the effective mass ratio (c) for a lattice depth Vb = 4.5-Er 
where Er = }i?k\/2m is the atomic recoil energy. The red 
solid line and blue dashed lines are for the first and second 
band case, respectively. Bottom: sketch of the density profile 
of a negative effective mass condensate in a trap potential 
U(x): the density is modulated by the optical lattice and 
pealcs at the maximum of U(x). " 



energy £'„(go). For slightly deep lattice (still far from the 
Mott insulator transition) it is sufficient to consider the 
first band only. This is the situation that we consider in 
the remainder of this paper. 

The validity of the negative effective mass description 
relics on the existence of a narrow momentum distribu- 
tion of the condensate relative to the central wave vector 
zbfci. This can be achieved by applying an initial ve- 
locity to the condensate, or alternatively by considering 
a condensate initially at rest and adiabatically switch- 
ing on a moving optical lattice realized by two counter- 
propagating fields of frequency difference 5ui = 2hkj^/m. 
This permits to prepare the condensate with a specific 
quasi-momentum [15[ such that the Hamiltonian (|10p de- 
scribes a quantum field of negative effective mass parti- 
cles trapped in the potential U{x) + £i{kL). 

The acceleration of particles with negative mass is op- 
posite to the direction of the forces to which they are 
subjected, so that stability occurs at the maximum of a 
potential U{x). We consider then the situation where a 
condensate of negative effective mass m\ j, is trapped in 
a potential U{x) that we approximate as an inverted har- 
monic potential of frequency il = ^J\U " {xi)£i {kL)\ / h 
with xi the position of its maximum. We expand the 
envelope field operator Ai_kL i^) on the basis of its eigen- 
functions £,e{x), with cigcnenergies huji = —hVl{l + 2) < 
0, as 



Ai.kLi^) = ^^ii{x)bi> 
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where the bosonic operator bi annihilates a condensed 



atom from the mode £. 

In contrast with the situation for positive masses the 
"ground mode" S,o{x), which has the largest population, 
has now the highest energy, But since the relative prob- 
ability for a particle to occupy state £ is proportional to 
the Boltzmann factor Pi oc exp[— fiw^/fesT], we conclude 
that the effective temperature T of a condensate with 
negative effective mass is also negative. We note that 
this situation is closely related to the scheme proposed in 
Ref. 18[ and recently demonstrated by Braun et al. [19| 
to achieve negative temperatures. In that case a tight 
optical lattice was used to establish the density profile of 
the ultracold atoms initially trapped in a loose harmonic 
potential U{x). That potential was then reversed from 
U{x) to ~U{x), following which the low-energy states of 
U{x) corresponded to the high-energy states of —U{x). 
In our approach, in contrast, U{x) is fixed and we move 
a shallow optical lattice to change the effective mass of 
the condensate from positive to negative. 

From that point on, the realization of a negative mass 
optomechanical system follows closely the approach pio- 
neered in Refs. [9|,ll0|, with an optical cavity field of wave 
vector kc <C k^ perturbing the center-of-mass motion of 
the condensate, but with important differences that we 
discuss in the following. 

For short times the depletion of the ground mode (,q{x) 
remains small and wc can describe it classically, bo ~ VN 
with TV the total atom number. The cavity-condensate 
coupling takes then the form of a multimode cavity op- 
tomechanical interaction 



be |20[ and hence a reversed optical output spectrum 
nc(— w) for the cavity field 
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where the effective optomechanical coupling coefficients 
are 



Gi = VNV / Co(a;) coa^ {kcX - e)U{x)dx. (14) 

Here T) is the single-photon potential depth of the cavity 
field of frequency Wc (including the frequency shift due 
to the condensate mean field) and annihilation operator 
c , ujmj == Wf — wq < are the oscillation frequencies of 
the "effective condensate mirrors," [9| and is a phase 
that depends on the position of the maximum of U{x) 
relative to the optical potential of the cavity field. We 
found numerically that coupling to the first excited mode 
^ = 1 can be made dominant by an appropriate choice 
of that phase, the remaining modes acting as a mechan- 
ical reservoir for that mode [11|. (Other situations with 
another single dominant mode (., or with two or more 
modes coupled with comparable strengths, can also be 
arranged.) 

For the inverted harmonic trap considered here the 
bath temperature is negative, as already mentioned. One 
can thus expect a reversed asymmetric displacement 
spectrum Sx{—iij) for the negative frequency oscillator 



This means that the 
optomechanical properties of a negative effective mass 
oscillator are reversed form the usual case. For example 
optical cooling is realized by a blue-detuned rather than 
red-detuned driving field, and stationary bipartite entan- 
glement is optimized in the red-detuned case [2l|, with 
the situation possibly even more interesting in the mul- 
timode case. When mechanical damping is weaker than 
the cavity decay, these negative frequency oscillators can 
serve as a negative temperature bath for a cavity field 
of positive frequency uic, so that it will exhibit gain as 
in usual laser theory - where a negative temperature is 
provided by the inversion of the active medium j22| . 

We are now in a position to discuss possible BEG op- 
tomechanical realizations of QMFS setups [5|. There 
are several ways to produce a pair of optomechanical 
oscillators with opposite effective masses. For example 
one could use a two-component condensate with only 
one component sensitive to the optical lattice, or ex- 
cite the condensate to the edges of two adjacent Bloch 
bands. Two-component condensates allow for the sensing 
and measurement of fields that couple differently to the 
two components, such as e.g. magnetic fields coupled to 
a spinor condensate or gravitational waves propagating 
along the system. 
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FIG. 2. Top: Possible arrangement for back-action-free field 
and force detection. Two condensates A and B are trapped 
along the axis of an optical cavity by the potentials Ua and 
Ub , respectively, with a moving optical lattice that drives con- 
densate B and adjusts its efi^ective mass to a negative value. 
Bottom: Relationship diagram for the setup in the "bare" 
(left) and "composite" (right) representations. The displace- 
ment of the composite subsystem D results in a change in the 
phase of the cavity field C that could be measured by homo- 
dyne detection, but the measurement back-action only affects 
the composite oscillator E. 

Figure [5] illustrates the situation where two cigar- 
shaped condensates A and B are trapped along the x- 



direction by the potentials Ua and Ub, respectively. A 
moving optical lattice interacts with the condensate B 
only and imparts it with a negative effective mass. The 
cavity resonance is sensitive to the motion of both con- 
densates, which form an effective "two-mirror" cavity op- 
tomechanical system with one positive and one negative 
mass. A QMFS is realized provided that the trapping po- 
tentials Ua and Ub and the optical lattice depth Vq are 
such that the optomechanical parameters of the two con- 
densates are identical except for the sign of their masses, 
and hence optomechanical oscillator frequencies. 

We assume for simplicity that a single mode is excited 
by the external perturbation to be measured, so that for 
each condensate a single term is dominant in Eq. (|13p . 
The effect of the weak coupling to other optical and en- 
velope modes can be described by a standard fluctuation 
and dissipation approach. The Hamiltonian part of the 
system dynamics is then governed by 

H = hwcC^c + hLj„i{a^a - b^b) + Gc^c{b + b^ + a + a'^) 
+ fA{a + a^) + fB{b + b^), (15) 

where a and b are the mode annihilation operators for the 
condensates A and B of mechanical frequencies com and 
—ujm, respectively, and G is the optomechanical coupling 
strength, see Eq. (fTB]). The parameters Ja and Jb, which 
are time-dependent in general, describe the coupling of 
the external perturbation to the condensates. 

The Heisenberg equations of motion for the cavity field 
and the composite displacement and momentum opera- 
tors Q = (a -h at -h 6 -h St)/2 and n = (a - at - 6 -h b'<)/2i 
are 
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where Ac = cUp — ujc is the cavity-pump detuning. These 
equations describe the motion of a particle driven by the 
difference between the strength of the external perturba- 
tions at condensates A and B. This causes a frequency 
shift of the cavity field that can be detected by interfcr- 
onictry. The key point here is that since these operators 
form a QMFS, [QjII] ~ 0, this measurement doesn't in- 
troduce any back-action to that QMFS, hence and the 
measurement is not subject to the standard quantum 
limit. Complementary conclusions are easily obtained 
for the second QMFS, characterized by the operators 
$ = (a -h at - 6 - St)/2 and P = (a - fit -h S - b^)/2i. 

Further insight into the underlying physics of this 
back-action-free measurement scheme by considering the 
quantum state dynamics. We assume that the system is 
initially uncorrelated, with the cavity field is in a coher- 
ent state and the oscillators A and B both in the ground 
state, 



As a result of the optomechanical interaction ([T5|) the os- 
cillators A and B become entangled with the cavity field 
C. However, when expressing the state of the system in 
terms of the composite oscillators D and E described by 
the operators {Q,P} and {$,n}, respectively, we find 
that the system does not suffer three-body entanglement 
between the subsystems C, D and E, but only two-body 
entanglement. Specifically we find, except for an unim- 
portant constant phase factor. 
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where 
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(A„(t) = ^=-^(./a + /b + 2Gn) (1 - e-'"-*) , 
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(19) 

and Qs = (/b — fA)/^A[- Equation ([T5|) shows that in 
contrast to the composite oscillator D, which becomes en- 
tangled with the cavity mode C, the composite oscillator 
E remains uncorrelated to the rest of the system. Rather, 
it evolves into a time-dependent coherent sate \ip{t))E 
that is independent of both the state of the optical field 
and the composite oscillator D. Its complex amplitude 
is fully determined by the steady displacement (|16p and 
the frequency ujm- 

The states (j)n{t))D are n-dependent coherent states of 
the composite oscillator D. They are similar to those en- 
countered in single-mirror optomechanical situations [23| , 
except for the important difference that the phase factor 
dependence depends on the steady-state displacement q^ 
of the oscillator E. That dependence makes it is easy to 
read out Qs without measurement back-action. For ex- 
ample for ujMt = 2rmr, n integer, the state of the system 
reduces to \aexp[—4im7TGqs])c'8\G)D^\0)E, that is, the 
composite oscillator D and E return to vacuum state - as 
do the oscillator A and B - while the cavity field becomes 
a coherent state whose phase could be easily measured 
by homodyne detection. 

Summarizing, we have proposed and analyzed a BEC- 
bascd optomechanical system that permits to isolate a 
QMFS with potential applications in back-action free 
measurements of feeble forces and fields. Our discussion 
centered on the use of quantum-degenerate atomic sys- 
tems whose spatial extent is large to the period of the op- 
tical lattice. One might ask if it is also possible to achieve 
a QMFS in ultracold, but non-condensed atoms. In that 
case every atom is localized within a single lattice well, 
but the phonon mode associated with the center of mass 
of the sample can still be delocalized and can enter the 
quantum regime via cavity optomechanical cooling [12 |. 
However thermal damping is now of the order of 7 '--^ ujm 



and is expected to significantly affect the measurement 
precision of the scheme. Future work will discuss in detail 
the effects of dissipation in the proposed system and con- 
sider specific applications, including magnetometry and 
gravitational wave detection. 
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